skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Law, Chi-Yan"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract We present JWST-NIRCam narrowband, 4.05μm Brαimages of the Sgr C Hiiregion, located in the central molecular zone (CMZ) of the Galaxy. Unlike any Hiiregion in the solar vicinity, the Sgr C plasma is dominated by filamentary structure in both Brαand the radio continuum. Some bright filaments, which form a fractured arc with a radius of about 1.85 pc centered on the Sgr C star-forming molecular clump, likely trace ionization fronts. The brightest filaments form a “π-shaped” structure in the center of the Hiiregion. Fainter filaments radiate away from the surface of the Sgr C molecular cloud. The filaments are emitting optically thin free–free emission, as revealed by spectral index measurements from 1.28 GHz (MeerKAT) to 97 GHz (Atacama Large Millimeter/submillimeter Array). But, the negative in-band 1 to 2 GHz spectral index in the MeerKAT data alone reveals the presence of a nonthermal component across the entire Sgr C Hiiregion. We argue that the plasma flow in Sgr C is controlled by magnetic fields, which confine the plasma to ropelike filaments or sheets. This results in the measured nonthermal component of low-frequency radio emission plasma, as well as a plasmaβ(thermal pressure divided by magnetic pressure) below 1, even in the densest regions. We speculate that all mature Hiiregions in the CMZ, and galactic nuclei in general, evolve in a magnetically dominated, low plasmaβregime. 
    more » « less
    Free, publicly-accessible full text available April 2, 2026
  2. Abstract We present James Webb Space Telescope (JWST) Near Infrared Camera observations of the massive star-forming molecular cloud Sagittarius C (Sgr C) in the Central Molecular Zone (CMZ). In conjunction with ancillary mid-IR and far-IR data, we characterize the two most massive protostars in Sgr C via spectral energy distribution (SED) fitting, estimating that they each have current masses ofm*∼ 20Mand surrounding envelope masses of ∼100M. We report a census of lower-mass protostars in Sgr C via a search for infrared counterparts to millimeter continuum dust cores found with the Atacama Large Millimeter/submillimeter Array (ALMA). We identify 88 molecular hydrogen outflow knot candidates originating from outflows from protostars in Sgr C, the first such unambiguous detections in the infrared in the CMZ. About a quarter of these are associated with flows from the two massive protostars in Sgr C; these extend for over 1 pc and are associated with outflows detected in ALMA SiO line data. An additional ∼40 features likely trace shocks in outflows powered by lower-mass protostars throughout the cloud. We report the discovery of a new star-forming region hosting two prominent bow shocks and several other line-emitting features driven by at least two protostars. We infer that one of these is forming a high-mass star given an SED-derived mass ofm*∼ 9Mand associated massive (∼90M) millimeter core and water maser. Finally, we identify a population of miscellaneous molecular hydrogen objects that do not appear to be associated with protostellar outflows. 
    more » « less
    Free, publicly-accessible full text available April 1, 2026
  3. Context. Evidence that the chemical characteristics around low- and high-mass protostars are similar has been found: notably, a variety of carbon-chain species and complex organic molecules (COMs) form around both types. On the other hand, the chemical compositions around intermediate-mass (IM) protostars (2M<m*< 8M) have not been studied with large samples. In particular, it is unclear the extent to which carbon-chain species form around them. Aims. We aim to obtain the chemical compositions of a sample of IM protostars, focusing particularly on carbon-chain species. We also aim to derive the rotational temperatures of HC5N to confirm whether carbon-chain species are formed in the warm gas around these stars. Methods. We conducted Q-band (31.5–50 GHz) line survey observations toward 11 mainly IM protostars with the Yebes 40 m radio telescope. The target protostars were selected from a subsample of the source list of the SOFIA Massive Star Formation project. Assuming local thermodynamic equilibrium, we derived the column densities of the detected molecules and the rotational temperatures of HC5N and CH3OH. Results. Nine carbon-chain species (HC3N, HC5N, C3H, C4Hlinear-H2CCC,cyclic-C3H2, CCS, C3S, and CH3CCH), three COMs (CH3OH, CH3CHO, and CH3CN), H2CCO, HNCO, and four simple sulfur-bearing species (13CS, C34S, HCS+, and H2CS) are detected. The rotational temperatures of HC5N are derived to be ~20–30 K in three IM protostars (Cepheus E, HH288, and IRAS 20293+3952). The rotational temperatures of CH3OH are derived in five IM sources and found to be similar to those of HC5N. Conclusions. The rotational temperatures of HC5N around the three IM protostars are very similar to those around low- and high-mass protostars. These results indicate that carbon-chain molecules are formed in lukewarm gas (~20–30 K) around IM protostars via the warm carbon-chain chemistry process. Thus, carbon-chain formation occurs ubiquitously in the warm gas around protostars across a wide range of stellar masses. Carbon-chain molecules and COMs coexist around most of the target IM protostars, which is similar to the situation for low- and high-mass protostars. In summary, the chemical characteristics around protostars are the same in the low-, intermediate- and high-mass regimes. 
    more » « less
    Free, publicly-accessible full text available December 1, 2025
  4. Abstract We present ∼8–40μm SOFIA-FORCAST images of seven regions of “clustered” star formation as part of the SOFIA Massive Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the turbulent core accretion (TCA) theory to derive key protostellar properties, including initial core mass,Mc, clump environment mass surface density, Σcl, and current protostellar mass,m*. We also carry out empirical graybody (GB) estimation of Σcl, which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python packagesedcreator, which is designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values ofMc∼ 30–200M, Σcl,SED∼ 0.1–3 g cm−2, andm*∼ 4–50M. The GB fitting yields smaller values of Σcl,GB≲ 1 g cm−2. From these results, we do not find evidence for a critical Σclneeded to form massive (≳8M) stars. However, we do find tentative evidence for a dearth of the most massive (m*≳ 30M) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function. 
    more » « less
    Free, publicly-accessible full text available June 3, 2026
  5. Abstract Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ 214μm observations of polarized thermal dust emission and high-resolution GBT-Argus C18O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion ofB-field orientations, we produce a map of the B-field strength of the IRDC, which exhibits values between ∼0.03 and 1 mG based on a refined Davis–Chandrasekhar–Fermi method proposed by Skalidis & Tassis. Comparing to a map of inferred density, the IRDC exhibits aB–nrelation with a power-law index of 0.51 ± 0.02, which is consistent with a scenario of magnetically regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions, and a dominant role forB-fields in regulating collapse. A clear alignment of magnetic field morphology with the direction of the steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with theB-field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation, and therefore they need to be accounted for in theoretical models of these processes. 
    more » « less
  6. Abstract We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from Atacama Large Millimeter/submillimeter Array 1.3 mm observations with a resolution of 0.″2 (∼1000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., H2CO, CH3OH, CH3OCH3), sulfur-bearing (SO2, H2S), and nitrogen-bearing (e.g., HNCO, NH2CHO, C2H3CN, C2H5CN) molecules. We discuss their spatial distributions, physical conditions, correlation between different species, and possible chemical origins. In the central region near the protostar, we identify three hot molecular cores (HMCs). HMC1 is part of a millimeter continuum ring-like structure, is closest in projection to the protostar, has the highest temperature of ∼300 K, and shows the most line-rich spectra. HMC2 is on the other side of the ring, has a temperature of ∼250 K, and is of intermediate chemical complexity. HMC3 is further away, ∼3000 au in projection, cooler (∼70 K), and is the least line-rich. The three HMCs have similar mass surface densities (∼10 g cm−2), number densities (nH∼ 109cm−3), and masses of a few solar masses. The total gas mass in the cores and in the region out to 3000 au is ∼25M, which is comparable to that of the central protostar. Based on spatial distributions of peak line intensities as a function of excitation energy, we infer that the HMCs are externally heated by the protostar. We estimate column densities and abundances of the detected species and discuss the implications for hot core astrochemistry. 
    more » « less
  7. Abstract We adopt the deep learning methodcasi-3d(convolutional approach to structure identification-3D) to infer the orientation of magnetic fields in sub-/trans-Alfvénic turbulent clouds from molecular line emission. We carry out magnetohydrodynamic simulations with different magnetic field strengths and use these to generate synthetic observations. We apply the 3D radiation transfer coderadmc-3dto model12CO and13CO (J = 1−0) line emission from the simulated clouds and then train acasi-3dmodel on these line emission data cubes to predict magnetic field morphology at the pixel level. The trainedcasi-3dmodel is able to infer magnetic field directions with a low error (≲10° for sub-Alfvénic samples and ≲30° for trans-Alfvénic samples). We further test the performance ofcasi-3don a real sub-/trans- Alfvénic region in Taurus. Thecasi-3dprediction is consistent with the magnetic field direction inferred from Planck dust polarization measurements. We use our developed methods to produce a new magnetic field map of Taurus that has a three times higher angular resolution than the Planck map. 
    more » « less
  8. Abstract We report high-resolution 1.3 mm continuum and molecular line observations of the massive protostar G28.20-0.05 with Atacama Large Millimeter/submillimeter Array. The continuum image reveals a ring-like structure with 2000 au radius, similar to morphology seen in archival 1.3 cm Very Large Array observations. Based on its spectral index and associated H30αemission, this structure mainly traces ionized gas. However, there is evidence for ∼30Mof dusty gas near the main millimeter continuum peak on one side of the ring, as well as in adjacent regions within 3000 au. A virial analysis on scales of ∼2000 au from hot core line emission yields a dynamical mass of ∼80M. A strong velocity gradient in the H30αemission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared spectral energy distribution (SED) analysis indicates a current protostellar mass ofm*∼ 40Mforming from a core with initial massMc∼ 300Min a clump with mass surface density of Σcl∼ 0.8 g cm−2. Thus the SED and other properties of the system can be understood in the context of core accretion models. A structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above ∼1Mfound from ∼0.1 to 0.4 pc around the source. This implies that a massive star can form in relative isolation, and the dearth of other protostellar companions within the ∼1 pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster. 
    more » « less
  9. Abstract We present ∼10–40μm SOFIA-FORCAST images of 11isolatedprotostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37μm imaging. We develop an automated method to define source aperture size using the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs with radiative transfer models, developed within the framework of turbulent core accretion (TCA) theory, to estimate key protostellar properties. Here, we release the sedcreator python package that carries out these methods. The SEDs are generally well fitted by the TCA models, from which we infer initial core massesMcranging from 20–430M, clump mass surface densities Σcl∼ 0.3–1.7 g cm−2, and current protostellar massesm*∼ 3–50M. From a uniform analysis of the 40 sources in the full SOMA survey to date, we find that massive protostars form across a wide range of clump mass surface density environments, placing constraints on theories that predict a minimum threshold Σclfor massive star formation. However, the upper end of them*−Σcldistribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher Σclconditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an ∼40 yr baseline. 
    more » « less
  10. Abstract We have obtained sensitive dust continuum polarization observations at 850 μ m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope as part of the B -fields in STar-forming Region Observations (BISTRO) survey. These observations allow us to probe magnetic field ( B -field) at high spatial resolution (∼2000 au or ∼0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis–Chandrasekhar–Fermi method, we estimate the B -field strengths in K04166, K04169, and Miz-8b to be 38 ± 14, 44 ± 16, and 12 ± 5 μ G, respectively. These cores show distinct mean B -field orientations. The B -field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B -field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. The B -field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B -field and not well correlated with other axes. In contrast, Miz-8b exhibits a disordered B -field that shows no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B -field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B -field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux. 
    more » « less